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Abstract 

Some magnetic properties of quark matter and a micro- 
scopic origin of the strong magnetic field in compact stars 
are discussed; ferromagnetic order is discussed with the 
Fermi liquid theory and possible appearance of spin den- 
sity wave is suggested within the NJL model. Implications 
of these magnetic properties are briefly discussed for com- 
pact stars. 

INTRODUCTION AND MOTIVATION 

Nowadays there have been many works about the QCD 
phase diagram. Here, we are concentrated in magnetic 
properties of quark matter and their implications on com- 
pact star phenomena. Let's begin with a simple question: 
what is and where can we expect magnetism in QCD. At 
low densities we may expect the appearance of pion con- 
densation in hadronic matter, where the classical pion field 
develops , followed by the specific spin-isospin order of 
nucleons 1 1 1 . In quark matter we shall see a non-uniform 
phase (called dual chiral density wave (DCDW) phase), 
accompanying the restoration of chiral symmetry, where 
the pseudoscalar condensate (qi^r^q) ^ as well as the 
scalar condensate spatially oscillates 0. Accordingly the 
magnetization of quark matter also oscillates like spin den- 
sity wave (SDW) in condensed matter physics. Further- 
more we may expect a ferromagnetic phase at some density 
region 0. 

On the other hand such magnetic properties should have 
some implications on the compact star phenomena. In par- 
ticular it has been well known that there is a strong mag- 
netic field in compact stars. The origin of such strong 
field is not clear even now, and it has been a long-standing 
problem since the first discovery of pulsars in early seven- 
ties. The recent discovery of magnetars seems to revive 
the problem again [4|. Their magnetic field amounts to 
O(10 15 G) from the P — P diagram. At present many peo- 
ple believe the inheritance of the magnetic field from the 
progenitor main-sequence stars or dynamo scenario due to 
the electron current. We consider here a microscopic origin 
of magnetic field by examining a possibility of spontaneous 
spin polarization in quark matter. 

Next, we consider a possibility of a non-uniform state 
in the vicinity of the chiral transition. Recently there have 
done many works about the non-uniform state at moderate 
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densities. We shall see that this phase exhibits an interest- 
ing magnetic property like SDW. 



FERROMAGNETIC TRANSITION 

The first study about the ferromagnetism in quark mat- 
ter has been performed by using the Bloch idea about the 
ferromagnetism of itinerant electrons |]5] . The mechanism 
is rather simple due to the Pauli principle: consider elec- 
trons interacting with each other by the Coulomb interac- 
tion in the background of the uniformly distributed posi- 
tive charge to compensate the electron charge. Then the 
Fock exchange interaction gives a leading-order contribu- 
tion. Then the electron pair with the same spin can effec- 
tively avoid the Coulomb interaction to give an attractive 
contribution to the total energy due to the Pauli principle. 
As the counter effect such polarized electron system costs 
more kinetic energy increase. So when the former effect 
becomes larger than the latter one, we can expect sponta- 
neous spin polarization. A perturbative calculation of quark 
matter interacting with the one-gluon-exchange interaction 
shows the transition to ferromagnetic phase at the order of 
nuclear density. Applying this idea to a star with solar mass 
and radius of 10km, we can roughly estimate the magnetic 
field of O(10 13 ~ 17 G). Thus we can feel that quark matter 
inside the core region may give an origin of magnetic field 
in compact stars. 

Recently we have studied the magnetic susceptibility 
Xm to get more insight about the properties of ferromag- 
netic transition within the Fermi-liquid theory [0. xm is 
written in terms of the quasiparticle interactions, 



Xm 



= fgpiM q \ 2 N(T) 



V 2 ) l + N(T)f a 
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where gn = f k ,_ fep df2k/47rg£)(k:) is the effective gyro- 
magnetic ratio, N(T) the effective density of states around 
the Fermi surface and f a the spin dependent Landau- 
Migdal parameter J7] 0. The spin susceptibility can be 
easily evaluated by using the OGE interaction. Then we 
can see that the Landau-Migdal (LM) parameters involve 
infrared (IR) divergences in gauge theories (QCD/QED). 
So we must take into account the screening effect at least 
to obtain the meaningful results. We have done it by cal- 
culating the quark polarization operator by the hard-dense - 
loop (HDL) resummation. As the results, we can see that 
the Debye screening for the longitudinal gluons surely im- 
proves the IR behavior, while the transverse gluons only 
receive the dynamic screening due to the Landau damping. 



Zero temperature case 

There are still left the divergences in the LM parameters 
at T = 0, but they cancel each other to give a finite xm- Fi- 
nally magnetic susceptibility is given as a sum of the con- 
tributions of the bare interaction and the static screening 
effect, 
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with k = m D /2k F in terms of the Debye mass, m D 
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Thus the screening effect 



g 2 pkp/2ir 2 , and C/ = 
gives the g 4 In g 2 term. 

To demonstrate the screening effect, we show in Fig. 4.3 
the magnetic susceptibility. We assume a flavor-symmetric 
quark matter, p u = Pd = Ps = Pb/3, an d take the QCD 
coupling constant as a s = g 2 /4tt = 2.2 and the strange 
quark mass m s = 300MeV inferred from the MIT bag 
model. Note that the screening effect is qualitatively dif- 
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Figure 1: Magnetic susceptibility at T = 0. Screen- 
ing effects are shown in comparison with the simple OGE 
case: the solid curve shows the result with the simple OGE 
without screening, while the dashed and dash-dotted ones 
shows the screening effect with Nt = 1 (only s quark)and 
Nf = 2 + 1 (w, d, s quarks), respectively. 



ferent, depending on the number of flavor Nf. The Debye 
mass is given by all the flavors, 



Fermi surface. Since the longitudinal gluons are short 
ranged due to the Debye screening, their contributions are 
almost temperature independent. Thus the main contribu- 
tion to the temperature dependence comes from the trans- 
verse gluons. Careful considerations about the quasiparti- 
cle energy show that quark matter behaves like marginal 
Fermi liquid, where the Fermi velocity and the renormal- 
ization factor vanish at the Fermi surface [0. Such behav- 
ior is brought about by the transverse gluons. The magnetic 
susceptibility is then given as 

(XM /XPauli) 1 = (XM/XPauli) 1 

C f9 2 u F (2k F + k 2 F m 2 + m i ) 2 



+ 



&k% 



2E 2 F 



+ 



72 

2 rri'2 \ 



T A In 



(4) 



with up = vp/Ep, where we can see the T 2 In T term ap- 
pears as a novel non-Fermi-liquid effect, besides the usual 
T 2 term, It should be interesting to compare this term with 
other ones in specific heat or the superconducting gap en- 
ergy. Furthermore, such logarithmic behavior also resem- 
bles the one by the spin fluctuations or paramagnons. Fi- 
nally the phase diagram is presented in Fig. [2] where we 
can also asses the importance of the non-Fermi-liquid ef- 
fect. 
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Figure 2: Magnetic phase diagram in the density- 
temperature plane. The open (filled) circle indicates the 
Curie temperature at kp = 1.1(1.6) fm _1 while the 
squares show those without the T 2 In T term. 



flavors 



(3) 



so that the «Tii(2/k) term changes its sign for k = 
m 2 D /2kp > 2. Thus we can see the screening^avors spon- 
taneous magnetization in large Nf. 

Non-Fermi-liquid effect at finite temperature 

We consider the low temperature case, T/fi <C 1, but 
usual low-temperature expansion can not be applied, since 
the quasiparticles exhibits an anomalous behavior near the 



Finally we present a phase diagram in Fig. [2] where we 
can estimate the Curie temperature of several tens of MeV. 
We can also see how the non-Fermi-liquid effect works for 
the ferromagnetic transition. 

MAGNETISM AND CHIRAL SYMMETRY 

Recently there have been appeared many studies about 
the non-uniform states in QCD, stimulated by the devel- 
opment of the studies about the exact solutions in 1+1 di- 
mensional models iflOl . The formation of the non-uniform 



states in quark matter have been studied in relation to chi- 
ral transition [11|. The appearance of density waves or 
crystalline structures has been an interesting possibility at 
moderate densities. Note that the appearance of the non- 
uniform phase is not special, but rather familiar in con- 
densed matter physics. In some cases it may exhibit an in- 
teresting magnetic property; the spin density wave (SDW) 
discussed by Overhauser is a typical example. In the previ- 
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Figure 3: Sketch of DCDW, where pseudoscalar density as 
well as scalar density oscillates along z direction. 

ous paper [2 1 we have discussed the possibility of a density 
wave, where pseudoscalar density as well as scalar density 
oscillate in harmony along one direction, which is called 
dual chiral density wave (DCDW). 
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Acos(0(r)), 
Asin(0(r)). 



(5) 



The chiral angle 0(r) is taken in the one dimensional form, 
0(r) = q • r. The amplitude A generates the dynamical 
mass M, M = — 2GA, while produces the axial-vector 
field for quarks, T3757 • V0/2 = T3757 • q/2. The single- 
particle (positive) energy is then given by 



E± = [El + q* /A±qy/pl + M*] l l\ 



(6) 



with E„ 



^Jp 2 + M 2 , depending on the spin degree of 
freedom. Accordingly the Fermi sea is split into two de- 
formed ones: one is deformed in the prolate shape and the 
other in the oblate shape. 

DCDW enjoys many interesting features. First, the sym- 
metry breaking pattern is Tp x Uqz(1) — > Up + Q3, which 
may be 1+1 dimensional analog of Skyrmion. Then the 
Nambu-Goldstone boson ("phason") has a hybrid nature of 
"pion" and "phonon". Secondly, a direct evaluation of the 
magnetization gives (012) oc cos(q • r), which means a 
kind of SDW. In this case quark matter can be regarded as 
a kind of liquid crystal endowed with two-dimensional fer- 
romagnetic order and one dimensional anti-ferromagnetic 
order. Note that magnetic field is globally vanished in this 
phase, but locally very strong. 

"Nesting " mechanism 

Here we discuss the mechanism for the formation of 
DCDW. There seems to be some confusions about it. In 
the references HI 21 authors emphasized the nesting effect 
(or Overhauser effect) for the essential mechanism of chiral 
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Figure 4: Energy spectra for p± = for q/2 > M . Solid 
(magenta) curves show the one for massive quarks, while 
dashed (blue) curves formassless ones, \p z ±q/2\. 



density waves, but there is little discussion about DCDW or 
other inhomogeneous phases. For 1+1 dimensional case, 
we can immediately see that q is given as q = 2[i for given 
chemical potential fi ifTUl . This is simply because the en- 



M 2 ± q/2 



ergy spectrum © is reduced to E^ 

and q is decoupled fromp^. Recall that the outstanding re- 
lation q = 2pp is held in the usual density wave like CD W 
or SDW in the one dimensional system, due to the nest- 
ing effect of the Fermi surface. We can see that the simi- 
lar mechanism works in the case of DCDW, but in some- 
what different manner from the usual one. For the case, 
M > q/2, Ep are only shifted ±q/2 from the free particle 
energy, so that formation of DCDW depends on the interac- 
tion strength like in the Stoner model. However, numerical 
calculation shows this is not the case: q/2 > M is always 
held in the DCDW phase. In Fig. H] we sketch the energy 
levels of the single quark energy for the case, q/2 > M. 
Note that mass is generated by the interaction with DCDW 
in this case. So E^ can be regarded as a consequence of the 
switching of the interaction with DCDW between massless 
quarks with relative momentum q. For massless quarks, the 
two levels cross each other at p z = for any q. Once the 
interaction with DCDW is present, mass is generated and 
two levels avoid the crossing with the energy gap, 2M, at 
p z = (magenta curves in Fig.|4|i. So if we choose q = 2[i 
and fill the levels up to pp = /i, there is always the energy 
gain due to the interaction with DCDW. In the three dimen- 
sional case, the simple relation is no more held, but we can 
expect some reminiscence. Actually we can numerically 
check that the similar relation is held in the three dimen- 
sional case. In the vicinity of the critical end point we have 
seen that the chiral correlation function x(l) diverges at fi- 
nite q of q ~ 2pp, but the effective mass is almost vanished 
in this situation. 

Thus we can understand DCDW with q/2 > M in 
terms of the "nesting" effect of the FErmi surface, while 
q smoothly increases from zero for RKC. Their situation is 
very different from our case: the opposite relation, \q/2 < 



M\, is realized in their case, and the level diagram looks 
very different from Fig. (|4). 

Deformed DCDW 

Here we generalize the original DCDW by taking into 
account the symmetry breaking effect with the current 
quark mass m c oc m\. Using a variational method, we can 
show that the functional form of the chiral angle in DCDW 
is deformed, satisfying the sine-Gordon (SG) equation, and 
thereby the allowed region of DCDW is extended. The sta- 
ble solution is then given in terms of the Jacobian elliptic 
function with modulus k, 

6 = tt + 2a,m (ml z/k,k) , (7) 

with the effective pion mass in medium, m* . To recover 
the original DCDW in the chiral limit, we must require the 
following relation, 

qk = m*7r/K (8) 
with the complete elliptic integral of the first kind K. Note 
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Figure 5: Relation of the modulus k and the parameter q. 
k — >• in the chiral limit. 

that we can recover the SG equation in ref . [13] and m* — > 
TOtt in the 1+1 dimensional case. 

CONCLUDING REMARKS 

We have discussed two kind of magnetic properties of 
quark matter separately, but a unified or comprehensive 
description of magnetism should be desired. Furthermore 
when we are interested in magnetism at moderate densities, 
we must take into account some non-perturbative effects 
explicitly. It may be one way to use the effective models of 
QCD, e.g. NJL model, to this end. 

Ferromagnetic order may have a direct implication on 
magnetic evolution of compact stars, but SDW order should 
also have some implications. DCDW may catalyse the el- 
ementary processes by providing extra momentum; for ex- 
ample it allows the quark /3-decay process as neutrino emis- 
sion during the thermal evolution. The magnetization is 
globally vanished there, but its fluctuation, (M 2 ), becomes 
large. Accordingly the local strong magnetic field may in- 
duce new QED processes. 



For the present it needs more studies about the proper- 
ties of the non-uniform states and relations among them. 
In particular, the comparison of DCDW and the real kink 
crystal in ifTTI is important, since they are typical structures 
in QCD, reflecting the different symmetries, U(l) vs Z%. 
More studies are needed including the symmetry breaking 
effect, thermal effect or model dependence. 

The existing region of magnetism may be overlapped 
with color superconductivity.lt is then interesting to elu- 
cidate the mutual relation in quark matter. Some works 
have been already done [14|, but more studies are needed, 
including unconventional mechanism of pairing. 
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